Quantum spin chains and Majorana states in arrays of coupled qubits 
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Several designs of inter-qubit coupling are considered. It is shown that by a combination of Joseph- 
son and capacitive coupling one can realize qubit interactions of variable spin content. Qubit arrays 
are discussed as models of quantum spin chains. In particular, a qubit model of the ID quantum 
Ising spin chain is proposed. A realization of unpaired Majorana fermion states in this system is 
considered. It is shown that Majorana states are represented by spin flip excitations localized on 
the chain ends. Using unpaired Majorana states in qubit chains for decoherence protected quantum 
computing is discussed. 



Quantum computers are machines that store informa- 
tion on quantum variables and that process that informa- 
tion by making those variables interact in a way that pre- 
serves quantum coherence . Quantum superconduct- 
ing circuits of nanoscale size became available recently for 
experiments which investigate novel quantum phenom- 
ena and their application to quantum computing |^~^). 
This research has led to the observation of Rabi oscilla- 
tions |j] , and to a demonstration of the superposition of 
macroscopic quantum states, the so-called "Schrodinger 
cat" states ■ The superposition character of these 

states was revealed juj by doing microwave spectroscopy 
experiments on the two quantum levels of superconduct- 
ing circuits consisting of three Josephson junctions com- 
prising the so so-called persistent current qubit [^] . The 
two states of this qubit have persistent currents of about 
a microamp and correspond to the motion of millions of 
electrons. 

In this article we discuss circuits consisting of many 
coupled qubits. As a simplest example of such a cir- 
cuit we consider an array of identical qubits with nearest 
neighbor couplings. Our interest in such circuits is two- 
fold. Firstly, qubit arrays represent a system in which 
novel types of qubit couplings with different 'spin op- 
erator content' and tunable control over the coupling 
strength can be investigated. Secondly, such arrays can 
be used to explore novel quantum phenomena with pos- 
sible applications in quantum computing. 

Even the simplest of these systems discussed below, 
the so-called 'quantum Ising spin chain,' is rather rich. 
Several interesting many-body states can be realized by 
varying qubit parameters and their coupling strength. 
Elementary excitations in some of these states are quan- 
tum solitons theoretically described as pseudofermions. 
We consider the phase diagram and point out several in- 
teresting regions in which the problem is exactly solvable. 
Also, as discussed below, the so-called Majorana fermion 
states can be realized as midgap states localized at the 
ends of the quantum Ising chain. In a recent article |T^ ] 
Kitaev outlined the possibility of using such states for 
reducing environmental sensitivity of quantum comput- 
ers. We show that the coupling of Majorana states to 
external time varying field can be weaker than that of a 



single qubit. 

The persistent current qubit comprises a super- 
conducting loop interrupted by three Josephson junctions 
(marked in Figj^(a) by an 'x') with the junction capac- 
itance C\£. The values of the three Josephson junc- 
tions coupling constants J\ and J2 are chosen so that 
the Josephson part of the Hamiltonian alone defines a 
bistable system which at the value of external magnetic 
flux $ = 0.5<I>o can be either in the right-hand, or in the 
left-hand current state. In this qubit, by choosing pa- 
rameters Ji t 2 and C'1.2 appropriately, the barrier in the 
phase space separating the right and left current states 
can be made low enough, so that tunneling between two 
classical states will take place. Also, detuning of external 
field from the value 0.5$o produces a bias on the right 
and left states, making one of them lower in energy than 
the other. 





a) 



FIG. 1. a) Qubit schematic shown with phases ±9 on su- 
perconducting islands, b) Josephson energy of the system 
as a function of the phase 9. Solid line: unbiased system 
with <E> = O.S^o; dashed line: system with finite detuning 
8<j> = $ - 0.5$ - 

It was shown in Ref. H] that the phase space coordinate 
which describes transitions between the right and left 
states is the relative phase of the two superconducting 
islands [Fig.[l](a)]. The potential energy of the system as 
a function of has two minima, = ±#0, as shown in 
Fig. [l](b). These minima are symmetric at $ = 0.5$o and 
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asymmetric at $ ^ 0.5<I>o. Tunneling takes place through 
the barrier separating the two minima. The resulting 
qubit Hamiltonian is 



H 



-Aa x - ha z 



(1) 



where A is the tunneling amplitude for the barrier U(9), 
and h = 2I P {<& — 0.5<I>o) is proportional to the detuning 
$ — O.S'&Oj and where I p is the circulating current. The 
energy of the ground state and the first excited states are 
given by 



E T = T Vh 2 + A 2 (2) 

In the notation of Eq. (||) the right and left current states 
correspond to qubit 'spin' up and down. The energy 
levels repulsion manifest in Eq. @ was directly observed 
recently by a spectroscopic technique p^O . 

The goal of this paper is to discuss a possibility to use 
superconducting quantum qubits to study novel quantum 
phenomena, control of qubit systems, and novel quan- 
tum devices. The question of how to make inter-qubit 
couplings of required strength, sign, and spin content is 
central for constructing logical gates out of qubits. We 
consider several basic couplings, such as zz, xx, etc., 
in one-dimensional qubit arrays coupled inductively, ca- 
pacitively, or by additional Josephson junctions. One- 
dimensional geometry is relatively simple to manufac- 
ture, as well as it is simplest to work on theoretically 
in order to have theory predictions compared with ex- 
perimental results. It is of interest to use arrays of cou- 
pled qubits to study novel quantum phenomena which 
can not be simulated either with present materials or 
present classical computers. In particular, qubit arrays 
can be used as models of quantum spin chains, one of the 
most basic systems in many-body physics, in parameter 
regimes that are not accessible with magnetic materials. 



I. INTER-QUBIT COUPLINGS; QUANTUM 
ISING PROBLEM 

Designing couplings is essential for being able to assem- 
ble qubits into circuits capable of performing computa- 
tion. For that, inter-qubit couplings must satisfy several 
requirements. First, the coupling strength should not be 
much less than the Rabi frequency of an individual qubit, 
so that the states of different qubits entangle with each 
other sufficiently quickly. Second, the coupling should be 
much weaker than the Josephson plasma frequency which 
sets the scale for the energy gap between the two lowest 
energy qubit states and states with higher energy. Third, 
the design should allow for different kinds of coupling, 
such as purely transverse, purely longitudinal, or more 
general couplings, as required by the particular scheme 
of processing quantum bits. Fourth, it should be possi- 
ble to turn the couplings on and off on sufficiently short 
time scale. Below we outline our approach to handling 
the problem of couplings. 



Conceptually, the simplest kind of inter-qubit coupling 
is of the longitudinal <j\(J 7 ^ form, which can be easily re- 
alized by magnetic interaction of qubits. This coupling 
results from magnetic field flux that the current in one 
qubit induces in the other nearby qubit. The form of 
this coupling in terms of qubit 'spin' operators, since the 
right and left current states have a definite ^-component 
of qubit 'spin,' is indeed erf erf . The sign of this coupling is 
such that it favors antiparallel spin state of two coupled 
qubits. However, the numerical value of the inductive 
coupling estimated for conducting loops of a few microns 
in size turns out to be at least an order of magnitude 
smaller that typical values of the tunneling amplitude A 
and detuning parameter h. Also, the inductive coupling 
value is preset by the system layout, often the lack of 
tunability, as the coupling strength is usually preset by 
the system layout. 

Because of that we propose another realization of the 
erf erf coupling which has values tunable in a wider range 
than that of inductive coupling. The coupling is achieved 
by a Josephson junction with a sufficiently large J shared 
by adjacent qubits as shown in Fig. ^|. In the regime when 
J ^ Ji,2 the phase drop across the large J junction is 
much smaller than the overall phase change. As a result, 
this coupling is gentle enough not to perturb the individ- 
ual qubit dynamics. However, the coupling strength t is 
of the order of J\il J , and thus it can be easily made of 
the same order as the single qubit parameters A and h. 
The sign of this coupling is positive. 




FIG. 2. ID array of qubits coupled by shared Josephson 
junctions: a realization of the erf erf interaction. 

As an example, we consider an array of qubits coupled 
by shared Josephson junctions, as shown in Fig. ^. The 
Hamiltonian for this system is 



(3) 



with positive t. Eq. (|3|) describes one of the basic many- 
body systems: an antiferromagnetic spin 1/2 chain with 
exchange constant t in external magnetic field with com- 
ponents h and A. We note that, since the inter-qubit 
coupling is of the erf erf kind, the 'exchange' interaction 
in the spin chain (^) is highly anisotropic. This should be 
contrasted with nearly isotropic spin exchange coupling 
encountered in magnetic systems, with usually only weak 
anisotropy possible due to spin-orbital coupling. 

There are several reasons for considering the one- 
dimensional qubit array shown in Fig. ^. One reason is 
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that magnetic measurements on this system can provide 
a very direct and simple test of inter-qubit couplings. To 
illustrate this point, we consider the system (||) in the 
limit A <C h, t. In this case the problem is reduced to 
the classical ID Ising model. Since the exchange coupling 
t is positive, the ground state in the absence of external 
field is antiferromagnetic. Weak external field \h\ -C t 
does not affect the ground state, while at high field all 
spins align. Magnetization as a function of the field h in 
this case is 



for h> 2t 
for \h\ < It 
for h < -It 



(4) 



This should be compared with the magnetization step re- 
ported recently for an array of isolated qubits JlCf] . The 
result (Q) means that due to the qubit interaction (||) 
the magnetization step splits into two distinct steps with 
relative separation equal to At. In a real system the mag- 
netization steps (P can be smeared due to randomness 
in qubit parameters, finite value of A, temperature, etc. 
Therefore, a demonstration of a magnetization curve sim- 
ilar to would indicate high level of qubit reproducibil- 
ity and control over qubit couplings. 

Another interesting fact is that the Hamiltonian (||) is 
nothing but the quantum ID Ising problem. Thus the 
qubit array shown in Fig. || provides a physical model of 
this system. The quantum ID Ising problem is known to 
have nontrivial ordered ground states and elementary ex- 
citations. The phase diagram of this system, sketched in 
Fig.||, can be easily obtained by mapping the problem (||) 
onto the classical 2D Ising model with anisotropic cou- 
plings and associating the strong coupling axis with the 
time direction jl3| . For generic parameter values h/t and 
A/t the problem is not integrable. However, on several 
lines indicated on the phase diagram there exist exact 
solutions. Different kinds of spin ordering in the ground 
state can be realized in this system. At zero 'external 
field' h there is (i) the ordered state with finite (erf) at 
|A| < t, and (ii) the 'disordered' state with finite (<rf } at 
|A| > t. The excitation spectrum has a gap at low ener- 
gies which vanishes at the critical points A = ±t (marked 
A and B in Fig. [§. 

It is instructive to consider the spin Hamiltonian (|^) 
at the JW line h = 0. In this case the problem can be 
solved by the Jordan- Wigner transformation jl^.|lt| 



2a] Oi 



i, of= (oi+4) n* 



(5) 



The operators m defined by (||) describe fermions obey- 
ing standard anticommutation relations, {ai,aj} = 0, 

|aj,<Zj| = Sij. In terms of the fermions a,i, the Hamil- 



tonian 



reads 



n = J2 t ( a i- 4) ( a i+i + 4+i) - 2Ac 4 a * ( 6 ) 



The transformation (||) and the resulting fermion prob- 
lem (|^) are applicable to both finite and infinite spin 
chains. 

h/t 



! class 



JW 



B 



A/t 



FIG. 3. The parameter space of the quantum Ising prob- 
lem (^). On the line 7 c iass the classical ID Ising model is 
realized. On the line JW the problem is solved by the Jor- 
dan- Wigner transformation (|E]), and has fermion excitations. 
On this line the problem is equivalent to the classical 2D Ising 
problem in the absence of external magnetic field, with the 
points A = ±t corresponding to the 2D Ising critical point. 
On the lines Z the problem is described by Zamolodchikov 
theory of critical Ising problem in magnetic field. 

We shall consider here the case of an infinite chain, 
when the Hamiltonian can be conveniently rewritten in 
the plane wave basis: 



H 



t [a p a- p e 



-ip 



h.c.) —2(A + t cosp) a p a p (7) 



This problem is readily diagonalized by Bogoliubov 
transformation yielding the dispersion relation of the 
form 



e(p) = 2\A + te l P\, H = ^e(p) b% 



(8) 



with the excitation wavenumber in the Brillouin zone 
— 7r < p < tt. The excitation spectrum (0) has a gap 



at low energies of the width E g 



2min A±i. The 



spectrum becomes gapless at A = ±t (the critical points 
A and B in the phase digram). 

Elementary excitations in this case are noninteracting 
Bogoliubov fermions. Let us remark that in terms of the 
original spin 1/2 variables (^) an energy excitation can be 
viewed as consisting of independent solitons which prop- 
agate through the chain preserving their identity even 
when other excitations are present. This situation is fa- 
miliar in integrable systems, classical and quantum. It 
may be of interest to consider using quantum solitons in 
qubit arrays as means of quantum information transmis- 
sion between distant parts of the qubit network. These 
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solitons can in principle be generated and detected lo- 
cally on individual qubits located at the opposite ends of 
the system. Such solitons can be used as information car- 
rier in essentially the same way as solitons in nonlinear 
optical channels. 



II. MAJORANA FERMIONS IN A QUBIT 
ARRAY 

In the theory of the quantum Ising model it is known 
Jl7[ that a rather natural representation of the spin prob- 
lem (||) can be obtained using Majorana fermions defined 
in terms of the Jordan- Wigner fermions as 



4 = 



1 

71 



1 V2: 



(9) 



The operators cy 2 are hermitian, = c", and satisfy 
the anticommutation relations 



{c?4+4c°}= 1 -5 aa ,6 jj ,, (c«) 2 



(10) 



The commutation relations indicate that each Majo- 
rana fermion represents a 'one-half of a Jordan- Wigner 
fcrmion, so that the latter can be viewed as a pair of two 
Majorana fermions. It was proposed by Kitaev jl2j that 
unpaired Majorana fermion states are better protected 
from environmental noise than conventional fermions. 
Thus it is of interest to look for realizations of these 
states. 

Here we demonstrate that the qubit array shown in 
Fig. ^| provides under certain conditions a realization of 
unpaired Majorana fermions. The Hamiltonian ([}]) for a 
chain of length n can be rewritten in terms of Majorana 
fermions as 



n 



3=1 



(11) 



The problem (11) can be diagonalized by an orthogonal 
2n x 2n Bogoliubov transformation of the operators Cj' , 
such that it brings the matrix 



(12) 



to a diagonal form. The spectrum of this problem con- 
sists of two parts. Most of the states have energies above 
the gap -Egap found for the infinite qubit chain. In ad- 
dition, there may exist two midgap states with energies 
within the gap near e = 0. The midgap states are local- 
ized near the edges of the system. 

To illustrate this point, we consider a semi- infinite 
chain, and look for an eigenvector of the semi-infinite 
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matrix ( |l2| ) with the eigenvalue e 
the form 



/ 1 \ 


A 



A 2 



A 



0. This vector has 



(13) 



\..J 

where A = A/t. For the state (|l3| ) to be normalizablc, 
one must have |A| < 1, in which case the normalization 
factor is A = (1 - A 2 )~ 1/2 . The condition |A| < 1 indi- 
cates that the midgap states occur in the ordered Ising 
phase t > |A| (the AB interval in the phase diagram). 

Note that only one Majorana species contribute to the 
midgap state, and so this is indeed an unpaired Majorana 
state. For a finite but long chain the midgap state on 
one end will be approximately given by (|l3|), whereas 
the state on the other end of the chain will be defined 
in a similar way in terms of the other Majorana fermion 
species. For a finite chain the two midgap states will be 
coupled (T^| by a matrix element of the order ~ A™. 

Let us discuss the meaning of these states in terms of 
the spin 1/2 formulation (§). It is instructive to consider 
the limit |A| <C 1, when the ground state of the spin 
chain is close to that of a classical antiferromagnet. In 
this case the width of the bulk gap is close to 2t, which is 
consistent with the energy required to reverse one spin in 
the bulk in the presence of coupling t to the neighbors on 
both sides. An approximately twice smaller energy t is 
required to reverse a spin at the end of the chain, which 
gives excitation energy in the middle of the bulk gap. 
Thus, each of the unpaired Majorana states at the ends 
of the chain is nothing but a single reversed spin, at finite 
A < 1 surrounded by slightly tilted spins on neighboring 
inner sites. 



tunneling 



FIG. 4. Tunneling between two unpaired Majorana exci- 
tations realized as midgap states of the ID quantum Ising 
problem. 

When the system is excited into one of the two midgap 
states, tunneling between two degenerate states will take 
place due to coupling of the midgap states on opposite 
ends across the chain. This tunneling process involves 
simultaneous coherent flipping of the spins at the chain 
ends accompanied by rearrangement of the tilted spins 
near the ends. It is noteworthy that this spin flipping 
does not change the net spin, since the two spins at the 
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ends flip coherently and in antiphase. Hence the (spa- 
tially uniform) background magnetic field fluctuations on 
the two ends are cancelled out and do not affect the tun- 
neling dynamics. This is consistent with the conclusion of 
Kitaev that unpaired Majorana states have weaker sen- 
sitivity to environmental time varying fields. 

III. TRANSVERSE INTER-QUBIT COUPLINGS 

Here we consider several interesting realizations of 
transverse couplings. One possible kind of transverse 
coupling is illustrated in Fig. |5j. In this scheme, supercon- 
ducting islands of adjacent qubits are capacitively cou- 
pled in such a way that the coupling capacitance C is 
of the order of the Josephson junction capacitances C\^- 
The physical reason for the coupling shown in Fig. || to be 
transverse is the following. When qubits are in a superpo- 
sition of the right and left current states, their quantum 
dynamics is described by tunneling between two minima 
of the potential U(9) in Fig. [|. Each tunneling event pro- 
duces a charge impulse on superconducting islands. By 
making the qubits capacitively coupled, one enforces the 
correlation of time-dependent charge fluctuations on dif- 
ferent islands, which is equivalent to having correlation 
of tunneling events on coupled qubits. The Hamiltonian 
for this coupling, as we shall see, is written in terms of 
the transverse spin operators a ± = \{o x ± i(J v ). 




FIG. 5. ID array of capacitively coupled qubits: a realiza- 
tion of transverse inter-qubit coupling described by the Hamil- 
tonian (0). 

To derive the form of the qubit coupling Hamiltonian, 
we consider two qubits coupled capacitively according to 
the scheme in Fig. |5| We represent the collective dynam- 
ics of the two qubits by motion in the two-dimensional 
phase space parameterized by the qubits' phases 9\, 9 2 - 
Potential energy is given by the sum of Josephson contri- 
butions U(9%) + U(92). In the absence of qubit coupling, 
their kinetic energy is m(9\ + #§)/2 with effective mass 
m being determined by the Josephson junction capaci- 
tances Ci 2 — see Ref. |]. Since the kinetic energy is 
isotropic, all tunneling amplitudes between four minima 
#1.2 = ±<?o are equal to each other. In this case, as shown 
schematically in Fig. ^, tunneling of the two qubits is in- 
dependent, and thus the tunneling Hamiltonian is written 
as a sum o~f + erf . Now, the capacitive coupling shown 
in Fig. [| is described by adding to the kinetic energy a 
new term m*(#i + ^2) 2 /8, where m* = (h/2e) 2 C. This 



term makes kinetic energy anisotropic by increasing ef- 
fective mass for the motion in the (f , 1) direction. As a 
result, when m* > m, all tunneling amplitudes are sup- 
pressed with the exception of the correlated transitions 
I Tl) I IT) illustrated in Fig. [| This process is de- 
scribed by the Hamiltonian of the form afa^ + a^a^- 
In the more general situation with m* ~ m, both the cor- 
related and independent tunneling processes take place, 
and the resulting Hamiltonian of the array has the form 

00 

H = J2 t(4^+i + °i°t+i) - (A< + haf) , (14) 

i— — 00 

where A is the individual qubit tunneling amplitude, and 
t is the correlated tunneling amplitude. Note that t 3> A 
when m* 3> m. The first term of Eq. ([TJ) can also be 
written as (t/2)(afaf +1 +ofof +1 ), i.e., the Hamiltonian 
( [Ti] ) defines the so-called XY spin model in an external 
field. 




FIG. 6. Correlated tunneling in two capacitively coupled 
qubits. 

Magnetization curve can be computed at A <C t by 
using the Jordan- Wigner mapping of the spin 1/2 prob- 
lem to a free fermion problem. In this case, as we discuss 
below, the mapping is constructed somewhat differently 
from the one used in the quantum Ising model. The re- 
sult is: 

{1, for h > t 

f arcsin(/i/*), lox\h\<t (15) 
-1, for h < -t 

In the transverse coupling case, instead of two abrupt 
magnetization steps (||) found for the qubit array with 
the zz coupling, we have a continuous magnetization 
curve with characteristic square root singularities at h = 
±t. 

To illustrate flexibility of the systems with transverse 
coupling, we consider an array with capacitive coupling 
built according to the scheme displayed in Fig. [?]. Like 
in the array in Fig. ||, effective qubit 'spin' couplings 
are purely transverse in this case. The analysis of cou- 
plings goes along the same lines as for the array in Fig. |^. 
Two decoupled qubits can be described by separable dy- 
namics in the two-dimensional phase space (6q, 6*2). The 
coupling changes kinetic energy by adding to it the term 
m*(0i — 02) 2 /2, where m* = (H/2e) 2 C. Note that in this 
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case, due to different orientation of coupled qubits, their 
coupling depends on the difference of the phases Q\ — 62 , 
rather than on the sum as in the above example. As 
a result, effective mass is higher for the motion in the 
(1, —1) direction, and at m* 3> m all tunneling processes 
except I ft) <-> I [I) are suppressed. According to Fig. ||, 
this process is described by the Hamiltonian of the form 
a i a 2 + a i a 2- A s before, for the more general situation 
of m* ~ m there are both the correlated and indepen- 
dent tunneling processes, and thus the Hamiltonian of 
the array in Fig. |^ is 

00 

H= ^t<i + - (A<rf + hat) ■ (16) 

i— — 00 

The first term of Eq. ( |l6| ) can also be written as 
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FIG. 7. Another realization of an array with transverse 
couplings described by the Hamiltonian (161). 

The difference between the problems (|Ll|) and ( |l6[ ) is 
most clearly displayed in the Jordan- Wigner representa- 
tion 

4 = 4 II °i > a i ^ a Ii a h °l = 2a t*i - 1- (17) 

To simplify matters we shall limit ourselves to the strong 
coupling case, when m» ^ m and thus t > A in both 
( |j"4| ) and ([To]). Then the fermionic representation of the 
Hamiltonians ( fu| ) and (16) is given by 



(i) : H = ^ t(a+a i+ i + a^^i) - 2/ia+a, , (18) 

i— — 00 
00 

(ii) : H' = i(a+a+ hl + a i+ iaj) - 2ha\ a t . (19) 



The spectrum of elementary excitations is obtained by 
going to the Bloch plane waves representation (in the case 
(ii) accompanied by Bogoliubov transformation). The 
result reads 



(i): e(p) = 2(tcos(p)-h) , 
(ii) : e(p) = ±2\jt 2 sin 2 (p) + h 2 



(20) 
(21) 



It follows from the result ( p0| ) that the ground states 
of the two systems are described by collective entangled 



states of all qubits in the chain, represented by ideal 
Fermi gas in the case (i) and by Dirac vacuum in the 
case (ii). In the case (i) the excitation spectrum is gap- 
less for not too large \h\ < t, whereas in the case (ii) there 
is a finite gap in the excitation spectrum equal to 4|/i|. 



IV. NOVEL QUANTUM PHENOMENA 

As we discussed above, there are several attractive 
schemes of designing qubit couplings. Besides allowing 
for control over coupling strength and its spin operator 
content, the proposed schemes are compatible with the 
plans to achieve tunability of couplings in time. This 
can be attempted, for example, by employing switchable 
Josephson junctions similar to the gated InGaAs super- 
conducting junctions [ fl8f which have very short switching 
times of the order of a few microseconds. Magnetization 
measurements on these arrays compared with theoretical 
results can provide useful information on qubit couplings. 

The one-dimensional arrays discussed above provide 
novel realization of spin 1/2 chains, one of the basic 
many-body physics models. Quantum spin chains are 
well studied theoretically and are known to be quite rich 
in properties |l5| , |l(| . Depending on the kind of coupling 
they can display different kinds of ordering: ferromag- 
netic or antiferromagnetic (both with algebraic correla- 
tions), spin liquid, Luttinger liquid, etc. Elementary ex- 
citations (quasiparticles) in these systems can have frac- 
tional quantum numbers and fractional (anyon) statis- 
tics. Depending on couplings, the excitation spectrum 
can be gapless or gapped, which should be manifest in 
dynamical properties, i.e. in the character of excitation 
transport along the chain. 

Experimentally, the only realization of spin chains that 
has been available so far is in quasi one-dimensional mag- 
netic materials. That, however, corresponds to a rel- 
atively small domain in the coupling parameter space, 
because spin exchange is usually isotropic or nearly 
isotropic, as a result of the weakness of spin-orbital in- 
teraction in solids. Also, dynamical phenomena such as 
excitation transport are quite difficult to access in mag- 
netic systems because available experimental techniques 
are limited to low frequency magnetization measurement 
and polarized neutron scattering. Thus realizing spin 
chains in the qubit arrays looks attractive from both 
points of view. First, in the qubit array couplings can 
be tuned to the desired form by the methods outlined 
in the previous section. Second, the dynamics can be 
directly probed by standard electric measurements. 

One can study excitation transport in the qubit ar- 
rays. Experimentally, this can be achieved by exciting 
the chain at one end, and measuring response at the op- 
posite end. It will be interesting to manufacture and 
compare properties of gapless and gapped systems, as the 
character of excitation transport has to be very different 
in the two cases. Moreover, in the case of the ID quan- 
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turn Ising problem (||), the dependence of the excitation 
energy (|8|) on the qubit coupling t indicates that for the 
array shown in Fig. || both regimes of excitation trans- 
port are available simultaneously. By changing the qubit 
coupling, one can sweep through the range of parameters 
in which the gap in the excitation spectrum (|) will grad- 
ually close and than reopen. The problem of excitation 
transport is also quite interesting theoretically. The chal- 
lenge here is that, in contrast with other one-dimensional 
Fermi systems such electrons in quantum wires, generic 
Hamiltonians of spin chains do not conserve the number 
of Jordan- Wigner fermions. For example, the external 
field term —haf in the ID quantum Ising Hamiltonian 
(|J) has a nonlocal form in the Jordan- Wigner represen- 
tation. This may in principle change the properties of 
elementary excitations. 

Looking somewhat ahead of experimental develop- 
ments, it is quite interesting to imagine and explore 
the possibilities that will become available if excitations 
in qubit arrays can indeed propagate ballistically, as 
suggested by the dispersion relation (||). One can ask 
whether it is possible to design qubit analogs of electron- 
based quantum wires and use them to manipulate ex- 
citations in the same way as it is done in conventional 
nanoscaled semiconductor structures. Doing this should 
include developing a 'battery' (i.e., a DC spin current 
source), as well as techniques of measuring the current 
flux of spin excitations. This is of interest because spin 
excitations, being neutral, are practically decoupled from 
external electric and (to a lesser extent) magnetic fields. 
Thus it may be possible to achieve levels of coherence un- 
available in electron systems (but perhaps comparable to 
what is possible for photons in quantum optics devices). 

Another interesting direction is to use qubit arrays to 
design unpaired Majorana fermion states. As the above 
discussion reveals, these states can be realized in the 
quantum Ising chain, i.e., in a qubit array with zz cou- 
plings. On the grounds of what has been discussed above, 
one can expect that these states can be probed spec- 
troscopically in a way similar to individual qubit states 
fiof . Using this technique one can, in principle, deter- 
mine decoherence time of Rabi oscillations realized with 
two Majorana states and verify theoretical expectation 
that unpaired Majorana states are protected from deco- 
herence. 
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